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ABSTRACT 

Given the radial velocity (RV) detection of an unseen companion, it is often of interest to estimate 
the probability that the companion also transits the primary star. Typically, one assumes a uniform 
distribution for the cosine of the inclination angle i of the companion's orbit. This yields the familiar 
estimate for the prior transit probability of ~ i?*/a, given the primary radius i?* and orbital semimajor 
axis a, and assuming small companions and a circular orbit. However, the posterior transit probability 
depends not only on the prior probability distribution of i but also on the prior probability distribution 
of the companion mass M c , given a measurement of the product of the two (the minimum mass 
M c sini) from an RV signal. In general, the posterior can be larger or smaller than the prior transit 
probability. We derive analytic expressions for the posterior transit probability assuming a power-law 
form for the distribution of true masses, dT/dM c oc M", for integer values —3 < a < 3. We show that 
for low transit probabilities, these probabilities reduce to a constant multiplicative factor f a of the 
corresponding prior transit probability, where f a in general depends on a and an assumed upper limit 
on the true mass. The prior and posterior probabilities are equal for a = — 1. The posterior transit 
probability is ~ 1.5 times larger than the prior for a — —3 and is ~ A/tt times larger for a = —2, 
but is less than the prior for a > 0, and can be arbitrarily small for a > 1. We also calculate the 
posterior transit probability in different mass regimes for two physically-motivated mass distributions 
of companions around Sun-like stars. We find that for Jupiter-mass planets, the posterior transit 
probability is roughly equal to the prior probability, whereas the posterior is likely higher for Super- 
Earths and Neptunes (10M® — 3QM®) and Super- Jupiters (3Mj up — 10Mj up ), owing to the predicted 
steep rise in the mass function toward smaller masses in these regimes. We therefore suggest that 
companions with minimum masses in these regimes might be better-than-expected targets for transit 
follow-up, and we identify promising targets from RV-detected planets in the literature. Finally, we 
consider the uncertainty in the transit probability arising from uncertainties in the input parameters, 
and the effect of ignoring the dependence of the transit probability on the true semimajor axis on i. 



1. INTRODUCTION 

Transiting planets have become the primary resource 
for characterizing the detailed properties of exoplan- 
ets. When complemented by RV measurements of the 
planet's orbital eccentricity, period, argument of perias- 
tron, and velocity semiamplitude, as well as by measure- 
ments of the stellar host's mass and radius, a photometric 
transit allows for the measurement of the mass, radius, 
density, and surface gravity of the planet. With these ba- 
sic parameters in hand, a transiting planet system is then 
amenable to a profusion of follow-up observations, which 
can then yield an impressi ve array of p hysical properties 
of the planet and star (see Winn||2011 for a review) . 

Transiting planets are currently discovered via two dif- 
ferent methods. First, planets that are initially discov- 
ered by RV measurements are monitored photometrically 
during the predicted time of inferior conjunction, in order 
to detect the small fraction that transit. In fact, the first 
identified transi ting planet HP 209458b was discovered 
in th is manner (Charbonneau et al. 2000 Henry et al. 
20001. Second, photometric transit surveys synoptically 



monitor large numbers of stars to search for eclipse sig- 
nals consistent with planetary-sized companions. These 
candidate transiting planets are then subjected to a bat- 
tery of follow-up observations in order to eliminate false 
positives and ultimately confirm the planet by measuring 
its RV signal and, in doing so, its mass. 
Because it is relatively easy and inexpensive to moni- 



tor large numbers of stars photometrically with the pre- 
cision needed to detect planetary transit signals, pho- 
tometric surveys for transiting planets can readily over- 
come the intrinsic rarity of transiting systems and effi- 
ciently identify large samples of candidates. Indeed, the 
vast majority of confirmed transiting planets were dis- 
covered in dedicated ground-based transit surveys (e.g. , 
Udalski et al.|[2002| |Konacki et ah] [20031 jAlonso et al 



2004| |McCullough et al |2006[ |Bakos et al. j|2007| |Collier 
Cameron et al.|2007WWeldrake et al.|2008[ |AIsubai et al 



|2011| |Siverd et al.||2012| . In addition, Kepler "has dis" 
covered over 2, 700 transiting planet candidates (Batalha 



et al.|2013 Burke et al.|2013 ) — an amount greater than 



that ot all previously- known exoplanets, although most 
of these lack precision RV confirmation and so lack pre- 
cision mass estimates. However, the cadence, photomet- 
ric precision, and total number of observations needed 
to achieve robust detections of transit signals typically 
requires monitoring relatively small fractions of the sky 
at a time; as a result, photometric transit surveys typi- 
cally identify transiting planets orbiting relatively faint 
(V > 10) stars. 

On the other hand, because of the relatively demand- 
ing observational requirements for precision RV measure- 
ments — in particular, high-resolution spectra, broad 
wavelength coverage, and a high signal-to-noise ratio — 
RV surveys must, in general, target bright stars in se- 
ries with relatively large aperture telescopes. Because 
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of these constraints, only a few thousands of stars have 
been monitored with precision RV. When combined with 
the rarity of planetary companions in general and tran- 
siting planets in particular, only eight of the planets orig- 
inally d etected via RV have been subsequently shown to 
transit (|Barbieri et al.l2007||Winn et al.||2011a| |Pemory| 
et~aT1|20l l) IChar bonneau et al-POOOj [Henry et al WMUft 
Sato et al.f|2005MBonhls et al.| | 20121 |Bouchy et al.|2005[ 
Moutou et al.|2009||Fossey et al.|2009[[Gillon et al.|2007p " 
IN evert heless, because they all orbit bright stars, these 
systems are the most amenable to follow-up programs 
and so are extraordinarily valuable. Indeed, these sys- 
tems are some of best-characterized planets outside our 
own solar system. 

Fortunately, it is likely that the both the sample size 
and diversity of transiting planets originally discovered 
by precise RV surveys will expand considerably in the fu- 
ture via several avenues. Amongst the currently-known 
RV planets, it is expected that there exist a handful 
of long-period transiting planets that have not yet been 
photometrically identified. The Transit Ephemeride s Re- 
finement and Monitoring Survey ( Kane et al. 2010a[ ) aims 
to identify these systems by first refining the orbits (and 
the estimated time of conjunction) of the most promis- 
ing systems (i.e. those with the highest transit prob abil- 



ities) and then foll owing these up photometrically ( Sea- 
groves et al. 2003). Even more promising, samples ot 
RV-detected planets have been and will continue grow- 
ing significantly. The measurement precision of optical 
RV surveys has steadily im proved, allowing for the detec- 
tion of first Neptune-mas s ( Mc Arthur et al.|2004 Butler 
et al.||2004[ |Lovis et al.||2006|), then Super-Earth mas"s 



( |Rivera et aI.|2005||Udry et al. j200_7j |Mayor et al 12 009 ) , 
and most recently Earth-mass ( |Pumusque et al.||20T2| 
companions, all of which ar e in trinsically more numer- 



ous than ga s giant planets (Howard et al 
et al. 2011). In the near 



2010 Mayor 



uture, the developm ent of 



dedicated near-IR precision RV instruments (e.g., Bean 



eFaLl[20T0l |Plavchan et"aLl[20T3l |Rayner fc PRVS Team 



2007 jMahadcvan et al.|2012p will enable surveys of large, 



previously-inaccessible samples of low-mass stars, and 
will also have enhanced sensitivity to low-mass compan- 
ions. Finally, massive parallel RV surveys, such as the 
Multi-object APO Radi al Velocity Exoplanet Large-area 
Surve y (MARVELS) (|Ge et al.||2008l lEisenstein et al. 



20111, may allow for surveys or giant planets orbiting 



somewhat fainter but considerably more numerous stars 
It is likely that substantial photometric follow-up re- 
sources will be necessary to identify the transiting sys- 
tems from these anticipated large samples of RV-detected 
exoplanets. For example, the CHaracterizing ExoPlan- 
ets Satellite (CHEOPS), planned for l aunch in 2017, is a 
satellite dedicated just to this purpose ( Feldt et al.|2007 ). 
It is important that these necessarily-limited pholomet- 
ric follow-up resources be allocated optimally to maxi- 
mize the yield of transiting planets. 

The most important input to optimizing follow-up ef- 
forts is an estimate of the transit probability P tr . The 
transit probability is defined as the probability, given 
both the properties of the companion that are inferred 
from RV measurements and the properties of the host 
star that are inferred from auxilliary measurements, that 
the orbit of the companion is inclined such that it passes 
in front of the host star from our perspective at some 



point during its orbit. The simplest and most commonly 
used estimate of the transit probability is (Borucki & 
Summers||1984| |Sackett||1999[ ), 



R* 
a 



(1) 



where i?„ is the radius of the host star and a is the semi- 
major axis of the companion's orbit. This estimate makes 
a number of assumptions that are widely appreciated, in- 
cluding: a circular orbit for the companion, a companion 
radius R c that is much smaller than the stellar radius, 
and a uniform distribution for the cosine of the inclina- 
tion angle i of the orbit. The effect of the finite size of 
the planetary companion can be easily included by, e.g., 
replacing i?* — > i?* + R c for grazing transits. The effect 
of eccentric orbits on the transi t probability has been 
explored in a number of papers dSeagroves et al."jj2003[ 
Barnesj2007"| |Burke|2008| [Kane fc von Braun||2008| ). Ec- 
centricity can boost the transit probability relative to the 
naive estimate in Equation]!] particularly for some favor- 
able geometries. HP 17156 b, with its ecc entricity of 0.67 
( |Fischer et al.|2007| [Barbieri et al.|2007|), and HP 80606 
b, wit h an e ccentricity ot .927 ( pNaet et al.|2001 Moutou 



et al.||2009[ | Fossey et al.||2009[), ar e dram atic examples 
of these effects. Beatty fc Seager (2010) found that a 
prior constraint on the inclination ot the host star can 



also boost the transit probability (see also Sackett||1999 
Watson et al.|20l"0 ), assuming the axis of trie companion 
orbit is aligned with the sp in axis of the star. Finally, 
Kane & von Braun (20081 showed that constraints on 
secondary eclipses can also affect the transit probability. 

However, there are some additional assumptions inher- 
ent in Equation [T] that are less widely appreciated and 
that can also affect the transit probability; we explore 
several of these here. In particular, the posterior proba- 
bility distribution of the orbit inclination angle i, given 
the detection of an RV companion, depends not only on 
the prior distribution of i but also on the prior distri- 
bution of the true mass of the companion M c . While 
the prior distribution of i is well-known (and is simply 
uniform in cosi), the prior distribution of M c is gener- 
ally not well-constrained — at least, not in the regimes 
of interest for exoplanet surveys. Although the effect of 
the prior distribution of M c on the transit probability 
has been noted and estimated previousl y for a few spe- 



cific cases (e.g., Wisniewski et al.||2012[), to the best of 



our knowledge it has not been explored in any detail. 
Previous authors have noted the dependence of the pos 
terior probability distribution of i (a nd thus M c ) on the 
prior distributions of i and M c (e.g., Ho fc Turner||20lf 



Lopez fc Jenkins|2012|. The fact that the transit proba- 



bility also depends on both these priors follows trivially 
from this result. In addition, here we explore two ad- 
ditional effects that we also believe have not previously 
been discussed in detail. First, we estimate the effects of 
uncertainties on the RV observables and host star prop- 
erties on the estimated transit probability. Second, we 
note that the semimajor axis in Equation [T] is properly 
the true semimajor axis (rather than the minimum semi- 
major axis), which is not known from RV measurements 
alone, and indeed depends on the orbital inclination for 
a given set of observables. 
The plan for this paper is as follows. In Section [2j we 
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review the observables for RV detected companions, and 
reiterate how these can be used to derive the familiar 
a priori transit probability from RV data. In Section 
[3l we derive the a posteriori transit probability for an 
observed minimum mass M and true mass distribution 
dN/dM c . In Section El we explore the posterior transit 
probability for powerdaw mass distributions, discuss is- 
sues of convergence for certain power-law distributions, 
and show that the posterior transit probability can be 
expressed as a simply scaling of the prior probability for 
low transit probabilities. We calculate posterior transit 
probabilities for different regimes of two model compan- 
ion mass distributions in Section [5] In Section [6j we 
compare the number of known transiting RV-discovered 
transiting planets with the number expected based on 
the naive prior transit probability. We also compare 
the posterior and prior transit probabilities for known 
RV planets, and identify planets whose posterior transit 
probabilities are both high and substantially higher than 
the corresponding naive prior probabilities. In Section 
[7J we discuss the effect of propagating parameter uncer- 
tainties on the estimate of the transit probability, as well 
as the effect of assuming the minimum semimajor axis 
to estimate the transit probability. We summarize our 
results in Section [51 



2. PROBLEM SET-UP 

Consider the RV detection of a faint companion, such 
as a single-lined spectroscopic binary. A well-sampled, 
high signal-to-noise ratio RV detection yields estimates 
of the orbital period T, semi-amplitude K, eccentricity 
e, the argument of periastron of the host star ui, and 
the time at some reference point in the orbit (e.g. the 
periastron). The semi-amplitude K is related to T, e, 
the host star mass M*, the companion mass M c , and the 
orbital inclination i via 



K 



2irG 



1/3 



M r sin i 



(M* + M c ) 2 /3 



(1 



-1/2 



(2) 



As is well known, the only model-independent physical 
quantity one can infer about a single-lined spectroscopic 
binary purely from direct observables is the mass func- 
tion M, which is defined as 



M 



(M c sini) 3 
{M* + Mc) 2 



K 3 T 
2nG ( 



.2^3/2 



(3) 



Generally, one can also estimate the mass (and radius 
J?») of the host star through a variety of methods. With 
an estimate of M* one can use the measured value of M. 
to determine M c as a function of i by solving the resulting 
cubic equation for M c . One can then also estimate the 
semimajor axis of the orbit as a function of i, 



( G[M, + M c (i)}\ 1/3 



4tt 2 



rp2/3 



(4) 



If one has independent reason to believe that M c <C 
Ms,, , then the "minimum mass" of the companion can be 



estimated directly by, 
Mq = M c sin i 

(—) 
\2irGJ 



1/3 



■ K 



(l-e 2 ) 1/2 (M*+M c ) 2 / 3 



1/3 



(1 - e 2 ) 1/2 M 



2/3 



(5) 



An excellent approximation for the condition that the 
companion transits the host is that its projected sepa- 
ration from the center of the host star at the time of 
inferior conjunction is less than or equal to the sum of 



R* and the radius of the companion R c (Kipping 2008 
Winn| 2011 ). The orbital separation of the companion 



from the host at the time of inferior conjunction is 

a(l-e 2 ) 

r c = , 

1 + e sin uj 

leading to the condition for a transit, 

r c cosi < i?* + R c . 

Or, in terms of a limit on the inclination, 

i?* + R c \ / 1 + e sin uj 
~a ) [ 1 - e 2 



cosi < 



(6) 



(7) 



(8) 



It is straightforward to demonstrate that, for 
isotropically-distributed orbit normals, the probability 
density distribution of cosz is uniform. In the absence 
of any information, we expect the orbits of binary sys- 
tems to be isotropic, i.e. to have no preference for a given 
orientation. Therefore, the a priori probability density 
distribution of cos i is expected to be uniform, and the a 
priori transit probability is given by 



Ptrfi — 



R* + R C 



1 



e sin uj 



1 -e 2 



(9) 



In the case of a circular orbit and R c <C R* , we recover 
the familiar transit probability Pt r ,o — R*/o. 

However, once one makes a measurement of the com- 
panion's minimum mass^ Mq — i.e. the product of M c 
and sini — the transit probability then depends not 
only on the prior probability density distribution of i 
but also on the prior probability density distribution of 
M c . That the transit probability depends on the prior 
on M c can be seen intuitively using the following exam- 
ple: Assume that one had prior knowledge that objects 
of mass < M m in did not exist. Then, if one detected 
a companion with minimum mass M Q = 0.1M mini one 
would be certain that the true mass M c was at least ten 
times larger than the minimum mass, and therefore that 
sin i < 0.1, or cosi > 0.995. Provided that the prior 
transit probability Pt r ,o < 0.995, the true a posteriori 
transit probability would be zero, since one would be 
certain that object did not transit. 

Wc derive the general expression for the a posteriori 
transit probability in the next section. Before moving 
on, however, we note one subtlety that we ignored in the 
previous discussion. The semimajor axis a in equation 
[9] is the true semimajor axis, which depends on M c and 

1 More precisely, once makes a measurement of the mass func- 
tion. 
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thus on i. Therefore, in order to estimate the transit 
probability, one must either solve simultaneously for a 
and the critical inclination for a transit, or use an ap- 
proximation to a. We return to this point in Section 

3. A POSTERIORI TRANSIT PROBABILITY 

Bayes' Theorem gives the conditional probability den- 
sity of event A given event B as 



P{A\B) 



P(B\A)P(A) 
~~ P(P) ■ 



(10) 



For our case, A = cos i and B = Mq. Then the probabil- 
ity density of cos i given a measurement of Mq is 



P(cosi|M ) 



P(M 1 cos i)P(cosi) 
P(M ) 



(11) 



where P(Mq) = J P(M \ cos i)P (cos i)d cos i and the in- 
tegral is over the allowed range of cosi. If we consider 
some maximum true mass M c>max , then the integral 
is evaluated from cos i = to some minimum inclin a- 
tion angle i min such that cosi mm = yl~ (sini mi „) = 
[1 - (Mq/M^) 2 ] 1 / 2 . 

The prior probability density of cos i is uniform; we 
write this as 



dr 

dcos 



= P(cosi) = const. 



(12) 



Here, T defines the rate or distribution, so dr/dcosz is 
the prior probability density with respect to cos i. The 
probability density of Mq given cosi, P(M \ cosi), is just 
the prior probability density of M p times the Jacobian 
between M c and Mq, i.e. 



P(M 1 cos i) 



dr dM c 
dM c dM 



dr 1 

dM c sin i ' 



(13) 



where dr/dAf c = P(M C ) is the prior probability density 
of the true mass M c . Thus, 



P(cosi|A/ ) = 



dr i 

dM c sin i 



Jo 



dr i 

dM c sin i 



dcos? 



(14) 



Finally, the probability that the companion transits is 
the (cumulative) probability that cosi is less than the 
critical value for a transit, which is 



P tr (cosz < X) 



Jo 



x d r i 

dM c sin I 



dcosi 



Jo 



^r^^dcosi 

dM c sin i 



(15) 



Here, we have defined the ratio of sum of the radii to r c 
i.e., the prior transit probability, as 



X = 



P* + P c 



P* +P C 



1 + e sin u> 

1-e 2 



(16) 



4. POSTERIOR TRANSIT PROBABILITIES FOR 
POWER-LAW MASS DISTRIBUTIONS 

4.1. Functional Forms 
Assume that dT/dM c — CM®, where C is a nor- 
malization constant. Then oT/dM c = C(4^) Q = 



CMq(1 — cos 2 i) a l 2 and Equation (15) becomes 



Ptr = 



Jo 



'(1 



cos- 



i)- 



~ d cos i 
2 d cos ? 



(17) 



It is clear from this expression that, for a = — 1 (a uni- 
form distribution in logM c ), the posterior transit prob- 
ability is equal to the prior transit probability. The inte- 
grand in the denominator is evaluated up to cos i m i n to 
allow for an upper mass limit on the prior distribution 
of M c . This is strictly necessary for a > 1 because this 
integral diverges. For a < 1, the integral converges, and 
so for simplicity we will take cos i m in — 1 f° r these cases. 
We discuss the effect of the choice of i m in (o r M c _ max ) 
on the transit probability for a > 1 in Section |4~2} 

Since X and cos i m i„ are never greater than one, we 
can re-write Equation (17) as 



Ptr — 



*[ 2 Pi(i 



g+i . 

2 • 



§;*") 



>.Fi(h, 



g+l . 3 



; cos z i r , 



0' 



(18) 



where 2F\(a, b; c; d) is the Gaussian (or ordinary) hyper- 
geometric function: 



2 F 1 (a,b;c;d) = ^ 



fe=0 



(a) k (b) k d k 
(c) k k\ 



(19) 



where (x) n is the Pochhammer symbol or rising factorial, 



if k = 




(20) 



The transit probability in equation [TT] can be written 
using analytic functions for specific values of a. Table 
[T] gives these functions for integer values of a between 
-3 and 3, and the solid lines in the top panel of Figure [T] 
shows Ptr versus X for these a values. Our posterior dis- 
tr ibutions for cos t agre e with the M c posteriors given M 
of Ho & Turner (2011 ) save for a = 1; using their nota- 
tion tor the numerator of the companion-mass posterior, 
$(M , X, a), we find that $(M , X, 1) = log(V A 2 - 1 + 
X) instead of $(M n ,AM) = log (2VX 2 - 1 + X) (cf. 
Equation 19 of |Ho fc Turner|201l| . 

Note that when x is a negative integer, (x)k is zero 
if k > — x + 1. The hypergeometric function in Eq. 
|19| thus has —6 + 1 terms when b is a negative in- 
teger. Therefore, when a is an odd negative integer 
(i.e. a = —1,-3,-5,...), the transit probability can 
be written as an odd polynomial in X of order |a| with 
(— a+l)/2 terms. For a = —1, the hypergeometric func- 
tions in the numerator and denominator of Equation [18] 
are unity, and thus the transit probability is simply A 
— i.e. the prior transit probability. 

Smaller a corresponds to a prior increasingly weighted 
toward smaller values of M c . For a given Mq, M c — 
Mq/ sin i oc 1/sini, so larger sini values are increasingly 
preferred for smaller a. As a result, smaller cosz values 
are more likely, and the transit probability is higher, as 
illustrated in Figure [T] 



We can also express Equation (17) as a function of 
period with Equation Q. We assume that = Mq, 
P* = Pq, and P c <C P*. The top panel of Figure[2]shows 
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X 

Figure 1. Posterior transit probability distributions as a function of the prior transit probability X for power-law true mass distribut ions 
with indices a = — 3, — 2, — 1, 0, 1, 2, and 3. For a > 1, we set M c ,max/Mo = 10 or sin i m i„ = 0.1 since the denominator of Equation j!7| 
does not converge for these values of a. For a = —1, the posterior is the same as the prior distribution, so Ptr = X. The solid lines show 
the exact posterior, while the dashed lines show the first-order Taylor series expansion centered at X = 0. For small X, the posterior 
transit probability is approximately a constant multiple of the prior probability. The bottom panel shows the fractional deviation between 
the full and approximate probabilities. 



the posterior transit probabilities for the same power- 
law true mass distributions as a function of period. The 
dependence on the power-law true mass distribution can 
create a factor of ~ 10 2 difference in the posterior transit 
probability at a given period for power-law exponents in 
the range a = —3 to 3. This disparity decreases for very 
small orbital periods since Pt r — > 1 as T — > 0. 

4.2. Dependence on Maximum M c for a > 1 



As noted previously, the denominator of Equation ( 17 1 
does not converge for cos i m i n — 1 when a > 1. In Figure 
[3]we show the sensitivity of the posterior transit probabil- 
ity to cosi m i n for a = 1, 2, and 3. We normalized Pt r to 
its value assuming M c , max /Mo = l(sini m i„ = 10. Since 
M c /M = J- = (1 - cos 2 i)~ 1/2 , increasing M Cymax cor- 
responds to increasing cosi m i„, and therefore lowering 
the transit probability. For a = 1, P tr diverges loga- 
rithmically, and therefore the precise choice of M c ^ max 
does not strongly affect the implied transit probability, 
for reasonable values. 

On the other hand, P tr diverges more strongly for 
a > 1, and in particular depends sensitively on M c max , 
and can be arbitrarily small. Accurate estimates of the 
transit probability of companions with minimum masses 
in such steeply-rising portions of the companion mass 
function therefore require reasonably narrow constraints 



to place a limit on M c ^ max . For example, companions 
around very massive main-sequence stars can often be 
ruled out based on the lack of evidence of flux from the 
companion or a second set of spectral lines (e.g., Fleming 
eIaLl[2012| ). 



4.3. Taylor Approximations 

The simple linear behavior for small values of X exhib- 
ited in Figures [T] and [2] invites us to make a first-order 
Taylor approximation f or P tr centered at X = 0. 
numerator of Equation (|18|) can be written as 



The 



X 



iFi 



1 3 



X' 



X 



k=0 



X 



2k 



l+2fc fc! 



-X[1 + 0(X 2 )}. 



(21) 



This allows us to express the Taylor approximation as 

P tr = f a X + 0(X 3 ). (22) 

In other words, the posterior transit probability for 
power-law priors on the companion mass is equal to the 
prior transit probability times a constant scale factor f a , 
up to third order in X. The scale factor f a depends 
only on a and i m - ln , and is simply equal to the inverse 



G 




Period (days) 



Figure 2. The top panel shows the posterior transit probability distributions as a function of period (solid lines) and their corresponding 
first-order Taylor approximations (dashed lines) for integer values of a between -3 and 3. We assume a solar mass and radius. For a given 
period, the top panel indicates that the posterior transit probability can differ by nearly a factor of 10 2 for different a. The difference in 
transit probabilities for exceptionally short periods decreases rapidly with decreasing period, converging to 1 for T — > 0. The bottom panel 
shows the fractional deviation of the first-order Taylor approximations to the full analytical transit probabilities. 



of the denominator in the full expression for the transit 
probability, 



fa 



Mil" 



1 3 



-;cos i r 



(23) 



For a < 1, we can use Gauss' Theorem to express f a 
in terms of Gamma functions, 



fa = 



2r(i- f) 



(24) 



It can be shown that this expression diverges for a > 1. 
Furthermore, for a <C — 1, we can use Stirling's ap- 
proximation that nl ~ \phrnn n e~ n a nd the identity 
n! = r(n + 1) to approximate equation (24 1 as 



fa 



-2a 
ire 



1 



a/2 



(25) 



An alternate expression for f a in this regime that is some- 
what less accurate, but also somewhat simpler is, 



fa 



-2a 



(26) 



definition of the exponential function, e x 

x/n) 7 



lim„ 



(1 



Figure ]4] shows the exact value of f a for a < 0, as 
the 1 



well as the two approximations give n in Equations 25 
and [26] We find that Equation [25] is an excellent ap- 
proximation to f a , deviating by < 10% Jot a < —1.7 
and by < 1% for a < 



-4.4. Equation 26 



is a some- 
what poorer, but still accurate, approximation, deviat- 
ing by < 10% for a < —2.5. Interestingly, the boost 
factor does not converge for a — > — oo, and thus the pos- 
terior transit probability can be arbitrarily large even 
for small prior transit probabilities, if one allows for a 
mass function that rises arbitrarily steeply toward small 
masses. This implies that detected RV companions with 
minimum masses near a large jump in the mass function 
toward lower masses will have quite high transit proba- 
bilities. 

The exact scale factors for several values of a are given 
in Table [T] We find the transit probability is boosted by 
~ 30% for a = -2, and by 50% for a = -30 However, 
the transit probability is reduced by ~ 64% for a = 
(/o = 2/tt), and by larger amounts for a > 0. 

The linear approximations to Pt r are plotted as dashed 
lines in the top panels of Figures [l] and [2j and the bottom 



This can be derived from either Equat ion 25 or directly 
from the definition of f a in equation [23] py using the 



2 We note that /_ 5 = 15/8 ~ 1.9, /_ 7 = 35/16 ~ 2.2, /_ 9 = 
315/128 ~ 2.5. 
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Figure 3. Posterior transit probabilities for a = 1, 2, and 3 as a function of M c /Mq, normalized to the corresponding transit probability 
with M c /Mo = 10. Increasing M c /Mq corresponds to increasing cos i m in- Ptr diverges only logarithmically for cos i m i n approaching unity 
for a = 1, but diverges more rapidly for a > 1, as illustrated for the specific cases a = 2 and 3. 



Table 1 

Posterior Transit Probabilities and Scale Factors 



Max X for % Error 

^TT% <5% < 10% 



Power-law Index a 



PtAX) 



-3 
-2 
-1 

1 
2 



! : (Xyjl - X 2 + arcsin X) 
X 

— arcsin X 

aictanh A' 

arctanh(cos * m ^ T1 ) 
tan imi n 



3/2 0.17 0.38 

4/tt 0.24 0.52 
1 

2/tt 0.24 0.52 

[arctanh(cosi mi „)] _1 0.17 0.38 

tani mi7l 0.14 0.31 

2[arctanh(cos i m i n ) + cot imin esc imin] 1 0.12 0.27 



0.52 
0.71 

0.71 
0.53 
0.44 

0.38 



panels show the fractional error between the full analyti- 
cal expressions and our Taylor approximations. The bot- 
tom panels indicate that the approximations are quite 
accurate for all but the smallest period orbits. In or- 
der to provide a more quantitative estimate of how well 
the Taylor expansions approximate the true posterior ex- 
pressions, we provide in the right columns of Table [T] the 
values of prior transit probability X (or R*/a for circu- 
lar orbits) for which the Taylor approximations deviate 
from the analytical expressions for P tr by 1%, 5%, and 
10%. For a < 0, the approximate transit probability ex- 
pressions are accurate to within 10% for most values of 
X; the 10% accuracy threshold moves to lower X with 
increasing a. For the eight transiting planets originally 
discovered by RV that were known at the time of pub- 



lication, X ranges from ~ 0.015 for GJ 436 b to - 0.29 
for 55 Cnc e, so the approximate transit probabilities for 
these planets have < 10% error (in fact, the approxima- 
tions for all but 55 Cnc e are within a 5% error). If we 
can tolerate errors of one part in 10, we can confidently 
use these approximate posterior expressions for most sit- 
uations. 

5. POSTERIOR TRANSIT PROBABILITIES FOR 
MODEL PLANET MASS DISTRIBUTIONS 

While the results from the previous section hold for 
power-law distributions, the true mass distribution of 
companions to stars of any type is not so well-behaved. 
For example, the brown dwarf desert is a well-known fea- 
ture in the mass function of relatively short-period com- 




Figure 4. Top panel: The first-order posterior transit scaling fa ctor f a as a function of a for negative values of a. The solid line shows the 
exact expressio n ob tained from numerical integration of Equation |23| with cos i m i„ = 1. The s hort -dashed line shows the approximation to 
fa in Equation |25| whereas the long-dashed line shows the simpler approximation in Equation |26| Bottom panel: The fractional deviation 
of the approximations to f a from the exact expression as a function of a. 



panions to solar- type stars (Marcy & Butler 2000). It 
represents a local minimum m this mass tunction, with 
the frequency of companions rising toward the stellar 
regime on the high-mass side of the desert and toward 
the planetary regime on the low-mass side of the desert 
( Grether fc Lineweaver|[2006 ) . 

The biggest obstacle to providing robust a posteriori 
transit probabilities is that the mass distribution of plan- 
etary companions is poorly known. Radial velocity de- 
tections yield only the minimum mass of the compan- 
ion, whereas planets discovered from ground-based tran- 
sit sur veys have severe and difficult-tq-quantify selection 



biases flPont et alj|2005| [Gaudi et al.||2005| |Gaudi||2005| 

Fressin et al. 2007). While Kepler has dramatically in- 



creased the number of known transiting planets, we do 
not have mass measurements for most of these systems. 

We therefore adopt synthetic, but plausible, physically- 
motivated planetary companion mass distributions that 



are derived from detailed planet formation simulations. 
We emphasize that we use these simply to illustrate the 
effects of more complicated mass functions on the transit 
probability; we do not claim to provide firm or robust 
predictions for these probabilities based on these mass 
functions. 

We consider ma ss functions predi cted by the planet for- 
mation models of Ida & Lin ( 2008 1 and Mordasini et al. 
(2012). To each of these, we also add stellar companions 
using a simple model based on the compilatio n of multi- 



ple stellar systems of Raghavan et al. (2010). Since the 



aridity of the brown dwart desert is poorly constrained 
we also consider two extreme variants of each mass func- 
tion: one with brown dwarf companions with mass be- 
tween the deuterium-burning and the hydrogen-burning 
limit as predicted by the planet formation theories, and 
one without such companions. We assume a deut erium- 
burning limit of M > 13Mj up ( jSpiegel et al.|2011 ) and a 
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hydro gen-burning limit of M > 0.07M© (Chabrier et al. 
20001. We also assume all companions are on circular 
orbits. We discuss the determination of these mass dis- 
tributions and analysis of the resulting posterior transit 
probability distributions in detail below. 

5.1. Discussion of the True Mass Distributions 

We took true masses of sub-ste llar companions from 
the population synthes is models of Ida & Lin ( 2004 1 and 



Table 2 

Scale Factors for Model Mass Distributions with Brown Dwarf 
Companions 



Mordasini et al.| (120121) . Although the companion mass 
function in these models depends on the period, for sim- 
plicity we included all planetary companions, regardless 
of period. 

For stellar-mass compani ons, we adopted a di s tribu- 
tion based on the resu l ts of Grether fe Lineweaver ( 2006 ) 
and |Raghavan et al. ( 2010 ~. Both groups examined a 
volume-limited sample of companions around Sun-like 
stars out to 25 pc from the Hipparcos catalog. Out to 
25 pc, the mass ratio of the companions around Sun-like 
stars, M c /M«, is nearly flat over the linear mass ratio 
range 0.2-0.95 (Raghavan et al. 2010). Assuming that 



the host stars are all of solar mass, then the distribution 
of companion masses is approximately flat in M c over 
the range 0.2M Q — 0.95M Q . For simplicity, we set our 
mass distribution to be flat in M c between 0.08M© and 
O.95M . 

We then combined the planetary and stellar mass dis- 
tributions by normalizing the different mass regions, as- 
suming that the planetary-mass companions from the 
model distributions with lMj up < M c < 13Mj up are 
twice as com mon as stellar-mass compa nions around 
Sun-like stars ( Grether fe Lineweaver|20"06 ). These com- 
bined mass distributions contain companions predicted 
by the planet formation models with masses in the brown 
dwarf regime. We also consider mass distributions with 
a completely arid brown dwarf desert. For these, we 
remove companions with masses between 13Mj up and 
0.07M Q . Thus, we consider four different ma ss distri- 
butions: two with masses from Ida fc Lin] (|2004[ ) and two 
with masses from Mordasini et al.~( 2012[ ), with one of 
each having a completely dry brown dwarf desert. 

For each of these M c distributions, we created a dis- 
tribution of minimum masses Mq as follows: For each 
true mass, we drew several values of cosi from a uni- 
form distribution, transformed these to sini, and then 
multiplied each true mass by the sini values to ob- 
tain minimum masses. Figure [5] show the final dis- 
tributions of true and minimum masses, binned in in 
0.25 dex bins in log mass. We show the distribu- 
tions obtained from the two planet formation models, 
each with and without a completely dry brown dwarf 
desert. To guide the discussion, we have subdivided the 
mass function into several (somewhat arbitrary) regimes: 
sub-Earths (10~ 8 M© - 0.1M©), Earths/Super- Earths 
(0.1M© - 10M ffi ), Neptunes (10M© - 100M©), Jupiters 
(100M© - 10 3 M©), Super- Jupiters (10 3 M© - 13M Jup ), 
brown dwarfs (13Mj up — 0.07M©) and stellar compan- 
ions (0.07M© - 1M ). 

The distributions all appear qualitatively similar and 
have several common features worth noting. In the stel- 
lar regime, the mass function falls toward lower masses 
as a power-law with a slope a ~ (by design). There 
is a local minimum in the brown dwarf regime, with 



Regime 



Mass Range 



Ida & Li 



LirR MordasinH 



Earths/Supcr-Earths 


O.lAf© - 


- 10M e 


-1.5 


1.12 


-1.4 


1.14 


Neptunes 


10M e - 


100M e 


-1.4 


1.07 


-1.7 


1.19 


Jupiters 


100M ffi - 


10 3 M e 


-0.9 


0.99 


-1.1 


0.99 


Super- Jupiters 


3M Jup - 


13M Jup 


-2.2 


1.42 


-1.1 


1.12 


Brown Dwarfs 


13M Jup - 


0.07M Q 


-0.9 


0.61 


-1.5 


0.86 


Stars 


0.07Af o 


- 1M 


-0.2 


1.80 


-0.1 


1.79 



a Mass distribution using Ida & Lin companion masses. 
b Mass distribution using Mordasini companion masses. 

the mass function rising toward lower masses continu- 
ing into the planetary regime. In particular, there is 
a relatively sharp rise for Super-Jupiters, with the mass 
function roughly behaving as a power law with a ~ —2 to 
a ~ —2.5 in this regime. The mass function for Jupiters 
is essentially flat in log(M c ), i.e., a ~ — 1. For lower- 
mass planets, the mass function again begins to rise, with 
a r~j —1.5 from the Neptune through the Earth/Super- 
Earth regime. Finally, the planet formation models pre- 
dict a peak in the mass function at or below an Earth 
mass, with a fall off for companions less massive than 
this peak. 

5.2. Posterior Transit Probability Calculations and 
Results 

We can now use the companion Mq values in the afore- 
mentioned mass distributions to estimate the posterior 
transit probabilities for planets in bins of AI Q . To pro- 
ceed, we assume that the host star is a solar radius and 
that the companion orbits with a semimajor axis of 0.1 
AU. We then calculate the transit probability of each Mo 
bin by examining the fraction of companions in that bin 
that transit — i.e. the fraction for which cos i < Rq / (0.1 
AU) . We then normalized the transit probability to the a 
priori transit probability i?©/(0.1 AU), thereby provid- 
ing an estimate of the factor / by which the a posteriori 
transit probability is scaled relative to the prior tran- 
sit probability, in analogy to f a for the power-law mass 
function priors. 

We plot the scale factor against mass in Figure [6j 
For most mass ranges, the scale factor / for a given 
M c sin i differs little between the different mass distribu- 
tions, with the most pronounced discrepancy happening 
in the brown dwarf regime below the deuterium-burning 
limit. In Table [2j we list the average scale factor val- 
ues in each mass regime for the two mass distributions 
with brown dwarfs. There are several important points 
to note about the distribution of scale factors. 

• In the Jupiter mass regime (100M© — 10 3 M©), the 
slope of the mass distribution is reasonably well- 
approximated by a power-law with index a w -1, 
so f_i = 1 and the prior and posterior transit prob- 
abilities are very similar. Consequently, using the 
posterior probability in this mass regime, rather 
than the prior, would not substantially change tar- 
get selection criteria for transit follow-up observa- 
tions (assuming the prior on the mass function in 
this regime is realistic). 

• In the Super- Jupiter regime (3Mj up — 13Mj up ), the 



10 



log(MJM 



los(MJM s ) 



Earths/ 
Super-Earths 



Earths/ 
Super-Earths 



f 



log(M c /M Q ; 
log(MJM e ) 



log(M c /M Q ) 
]ag(MJM e ) 



Earths/ 
Super-Earths 



Earths/ 
Super-Earths 



I 



log(M c /M ) 



log(M c /M Q 



Figure 5. Mass distributions of companions around Sun-like stars using planetary companions from Ida & Lin (top) and Mordasini 
(bottom), with (left) and without (right) brown dwarf companions in the deuterium-burning to hydrogen-burning mass range. The 
true mass distributions are shown as filled histograms, whereas the minimum mass distributions are outlined. The vertical dashed lines 
divide the plots into the sub-Earth, Earth/Super-Earth, Neptune, Jupiter, Super- Jupiter, Brown Dwarf and stellar mass regimes as 
described in the text. Values of a for the various regimes are as follows, starting with Earths/Supcr-Earths and increasing in mass: 
a « -1.5,-1.4, -0.9, -2.2, -0.9, and - 0.2 for the top-left panel; a as -1.5, -1.4, -0.9, -2.5, -0.8, and - 0.1 for the top-right panel; 
a RJ —1.4, —1.7, —1.1, —1.1, —1.5, and —0.1 for the bottom-left panel; and a —1.4, —1.7, —1.1, —1.6, —1.3, and —0.1 for the bottom-right 
panel. 



transit probabilities are boosted by ~ 20 — 50%. 

• The posterior transit probabilities for Neptuncs 
and Super-Earths are also boosted relative to the 
prior transit probabilities. 

• The posterior transit probabilities for companions 
with minimum mass near the high-mass end of the 
brown dwarf desert have very low posterior transit 
probabilities and are thus poor targets for transit 
follow-up searches. These systems are simply much 
more likely to be low-mass stars seen at relatively 
low inclination angle i. 

Together, these results suggest that RV detected planets 
with M c sini in the mass range of Earths, Super-Earths, 
Neptunes, or Super- Jupiters are generally better-than- 
expected targets for transit follow-up, whereas compan- 
ions in the BD desert are poor targets for transit follow- 
up. 

6. APPLICATION TO KNOWN RV SYSTEMS: 
IDENTIFYING PROMISING SYSTEMS FOR 
PHOTOMETRIC FOLLOW-UP 

In this section, we apply the posterior transit prob- 
abilities derived from our model mass distributions to 



known RV planets in order to identify specific systems 
that might be significantly better candidates for transit 
follow-up observations than would be expected based on 
the naive prior transit probability. 

Before doing so, however, it is worth asking whether 
or not there is any evidence that the planet mass dis- 
tributions we have adopted have any correspondence to 
the true mass distribu tion of planetary companions (c.f. 



Mordasini et al. 



2009) 



Estimates of the true mass distribution of Jupiter-mass 
planets discovered by RV surveys generall y find distribu- 
tions that are roughly flat in log M c (e.g. Watson et al. 



20101, consistent with mass functions used here and sug- 



gestmg that the posterior and prior transit probabilities 
for such companions should be similar. It is wo rth noting 
that the first transiting planet, H P 209458 b (Charbon- 
neau et al.|2000 Henry et aL]|2000[ ), was found "right on 



time" based on adopting the prior transit probability for 
hot Jupiters of ~ 10%, and given the ~10 hot Jupiter 
systems that were known at the time. 

It is also interesting to note that there are two transit- 
ing Super- Jupiters and three transiting Jupiters among 
the planets orbiting solar-type stars that were originally 
discovered by RV, despite the fact that super-Jupiters are 
much rarer than Jupiter-mass companions among short- 
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Figure 6. Scale factors / (ratio of posterior to prior transit probability) as a function of minimum companion mass for the model 
mass distributions. The upper-left panel corresponds to the distributions with Ida & Lin planetary masses while the upper-right panel 
corresponds to Mordasini planetary masses. The bottom panels compare scale factors as a function of mass between the Ida & Lin and 
Mordasini distributions with (bottom-left) and without (bottom-right) brown dwarfs. 



period systems. This is consistent with the prediction 
of the mass functions adopted here that the posterior 
transit probability of Super- Jupiters is higher than the 
posterior transit probability of Jupiters. 

While this anecdotal evidence based on the rates 
of transiting Jupiters and Super-Jupiters is intriguing, 
there are of course other factors that strongly affect the 
transit probability, including the semimajor axis, argu- 
ment of periastron, and host radius. To account for these 
effects, we compared the number of RV-detected plan- 
ets that have been identified to transit to the expected 
number of transiting planets based on the naive prior 
transit probability. We took our exoplanet sam ple from 



the Exopla net Orbit Database on 19 April 2013 (Wright 
et al.|2011| ). We included only the 380 RV planets in the 



database whose eccentricity e, semimajor axis a, angle of 
periastron u> and host star's radius i?* were listed. We 
also included the eight planets first detected by RV and 
subsequently found to transit. For each RV planet, we 
found the prior transit probability from its orbital pa- 
rameters and using the general form of the prior transit 
probability given in Equation ([£]), assuming i?* ^> R c . 
We then partitioned the planets into three 1-dex bins 
of log(M c sinz). We estimated the expected number of 
transits in each bin by adding the prior transit proba- 
bilities of the planets in that bin. We divided the actual 
number of transiting planets in each bin by the predicted 
number of transiting planets. 



This ratio of the number of known transiting planets 
to the number predicted based on the prior transit prob- 
ability is shown in Figure [7] Note that the ratio plotted 
is strictly a lower limit to the actual ratio of transiting 
planets to predicted number of transiting planets, since 
there may exist transiting planets among the sample of 
RV detected systems that have not yet been identified 
to transit. In particular, identifying or definitively ex- 
cluding transits of long period systems is quite difficult 
because of the infrequent transit opportunities and gen- 
erally large uncertainties in the predicted times of infe- 
rior conjunction. With this in mind, we also show the 
ratio after restricting to RV systems with semiamajor 
axes a < 0.5 AU and a < 0.2 AU. Presumably, the sys- 
tems in these subsets have been more thoroughly vetted 
for transits, and thus the sample of transiting planets is 
more nearly complete. Because of the small number of 
known transiting systems and the likely incompleteness, 
we caution against drawing any strong conclusions. Nev- 
ertheless, for the bins in which there are a non-zero num- 
ber of known transiting systems, the results are gener- 
ally consistent with the expectations from the mass func- 
tions adopted here. In particular, the observed number 
of transiting planets with minimum mass between Mj up 
and 10Mj U p (M c s\ni = 10 25 M ffi - 10 35 M ffi ) transiting 
planets is apparently higher than one would expect based 
on the naive prior transit probability. 
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Figure 7. The ratio of transits to the expected number of transits as a function of M c sini, plotted for all planets in the sample (solid 
line), planets with semimajor axes a < 0.5AU (dashed line) and planets with semimajor axes a < 0.2AU (dotted line). 



6.1. Posterior-Prior Comparison and Transit 
Candidate Identification 

To compare the prior and posterior transit probabili- 
ties for known RV planets, in Figure [8] we plot the pos- 
terior probability versus T, following Kane & Von Braun 
(2008), as well as the ratio of the posterior and prior tran- 
sit probabilities. The prior transit probabilities were es- 
timated with Equation [9] assuming R p -C R* , the values 
of e, L), i?*, M*, and a as given by the Exoplanets Orbit 
Database]^] We determine posterior transit probability 
scalings using both the Ida & Lin and Mordasini model 
mass distributions that include brown dwarfs. To find 
the posterior transit probability scaling for a given RV 
system, we simply choose the value of the posterior tran- 
sit probability scaling for the 0.25 dex bin in log(M c sini) 
in which the RV companion is located in Figure [5j 

The overall trend of the transit probabilities versus 
T is similar for all cases, and simply reflects the scal- 
ing P tr oc X oc T~ 2 / 3 . To guide the eye, in Figure 
[8] we plot the simple prior transit probability for cir- 
cular orbits, a star with a solar density, and a planet 
with a small mass and radius compared to the star: 

P tr = (3TT/Gy/ 3 pQ 1/3 T- 2 / 3 . Deviations from this fidu- 
cial probability are due to eccentricity, orientation of the 
orbit (trj), variations in p*, and finally, in the case of the 

3 The semimajor axes in this database were derived using Ke- 
pler's Third Law and assu ming that the true companion mass is 
equal to its minimum mass ( |ButIer etT al. 2006 Wrigh tTet al,|20TT| ) 



posterior probabilities, variations in the posterior scal- 
ing factor. The vertical grouping of systems with large 
periods (10 2 days < T < 10 3 days) and relatively large 
transit probabilities are planets orbiting giant stars with 
large radii. Although the transit probabilities are large 
for these systems, any transits in such systems would be 
very shallow and very long. Specialized methods are gen- 
erally need ed to detect such t ransits; these are discussed 
in detail in |Assef et al.| ( |2009[) . 

The transit scale factor (ratio of posterior to prior tran- 
sit probabilities) inferred from the Ida & Lin model shows 
considerably more variance than that inferred from the 
Mordasini model. This simply reflects fact that the Mor- 
dasini mass distribution is smoother than that of the Ida 
& Lin distribution in the planetary regime (see Fig. [5]). 
Also, because the Mordasini distribution of companion 
masses increases nearly monotonically toward lower mass 
down to ~ M0 , the scale factors inferred from this model 
are never substantially less than unity for the known ex- 
oplanets. On the other hand, the Ida & Lin mass dis- 
tribution shows a clear local minimum near ~ 30 M ffi (a 
consequen ce of the "planet desert" discussed in |Ida fc| 
|Lin||2004| , and therefore the inferred scale factors for 
planets with minimum mass somewhat larger than this 
minimum are significantly less than unity. For the Ida & 
Lin model, we find that the posterior probability is less 
than the prior probability for ~ 19% the planets in the 
sample, as opposed to ~ 34% for the Mordasini distribu- 
tion. 

More interesting, however, is the fact that we infer a 
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posterior probability that is markedly higher for some 
planets, particularly for the Ida & Lin model. These are 
the ones we will focus on here. To identify such promising 
transit candidates, we begin by eliminating giant star 
hosts with R* > 2.5Rq (marked with plus signs in Figure 
[8 1 from consideration. From the remaining planets, we 
choose planets from either model that have both a high 
posterior transit probability of P tr > 0.1, and a high 
prior transit scale factor of / > 1.2, i.e. for which P tr > 
l.2P tr ,Q. These planets lie above the dotted lines in both 
the top and bottom panels of Figure [5J and are marked 
with an X. 

These cuts leave us with a sample of fourteen planets, 
which are listed in Table [3] Thirteen planets pass the 
cuts from the Ida & Lin distribution, whereas only HD 
47186 b passes from the Mordasini distribution. None 
of the candidates pass the cuts from both distributions. 
Among these fourteen, 55 Cnc e and HD 17156 b have 
already been shown to transit, whereas transits in HD 
40306 b, HIP 14810 b, and t Boo b have been con- 
clusively ruled out. The sum of the posterior transit 
probabilities of the remaining candidates is ~1.15, sug- 
gesting that one transiting system may lurk amongst 
these systems. Detecting or excluding transits from most 
of these systems will be challenging from the ground, 
as most have predicted depths of ^0.05 — 0.1%, based 
on the minimum mass and a mass/rad ius relation of 
R c = i? ffi (M c /M ffi ) - 53 ([Weiss et al.|2013|. 



7. 



ADDITIONAL ISSUES AFFECTING THE 
TRANSIT PROBABILITY 



We now turn our attention to two issues that, to our 
knowledge, have not been previously discussed and can, 
in principle, subtly affect transit probability estimates. 
First, we consider the uncertainty in the transit proba- 
bility arising from uncertainties in the input parameters 
that are used to estimate the transit probability. Second, 
we consider the effect of the dependence of the semima- 
jor axis (and so the transit probability) on the orbital 
inclination. In both instances, we find that the quan- 
titative changes to the estimated transit probability are 
likely to be negligible in most instances, although they 
can be larger and therefore more important in some spe- 
cial cases. 

We begin by writing down the expression for the poste- 
rior transit probability in the limit of small transit prob- 
abilities (i.e., X <C 1), in which case we can write Pt r as 
a constant factor / times the prior transit probability)^] 



P tr ~fX = f 



where we have defined 



R* + R c 



?(e,w). 



e sin u) 



1-e 2 



(27) 



(28) 



It is instructive to deconstruct this expression for the 
transit probability in terms of the quantities that can be 
measured directly from the RV data and those that must 
be inferred or assumed from external information. The 
relevant RV observables are T, e, and u>. Using Newton's 

4 Formally, we only demonstrated this to be true for power-law 
priors for M c , but we expect it to be generally true for any prior 



form of Kepler's Third law and defining the bulk density 
of the star p* = 2>M*/(AitR^) : the transit probality can 
be written as 



Ptr — 



1/3 



fp- 1/3 (l + q )-^(l + r)T 



1/3/ 



-2/3 



(e,u>) 



(29) 

where we have defined the mass ratio q = M c /M„ and the 
radius ratio r = i? c /i?*. Thus, in the limit that r -c 1 
and q <C 1 the only parameter of the star that enters into 
the transit probablity is the density p*. In fact, one can 
write, 



Ptr = (3 7 r) 1 / 3 /(l + g)- 1/3 (l + r) 



^dyn 



-2/3 



5(e,w), 
( 3 °) 

suggesting that the more fundamental physical quantity 
is, in fact, the free-fall or dynamical time of the star 

tdyn = (G/3*) -1 / 2 . 

It is useful to examin e an d classify the parameters 



that enter into Equation (29). The RV observables T, e, 
and oj will have uncertainties associated with them, and 
these uncertainties may be correlated. The stellar den- 
sity p* must be inferred from external information, i.e. 
from the stellar temperature, surface gravity, and metal- 
licity derived from a high-resolution spectrum (perhaps 
combined with a parallax measurement) or directly from 
astroseismology. Both the RV parameters and stellar 
densities will have measurement uncertainties associated 
with the m. W e discuss the effect of these uncertainties 
in Section [TT Finally, the posterior scaling factor / and 



the radius ratio r must simply be assumed. These may 
be uncertain, but these are not statistical uncertainties in 
the traditional sense; rather, they are systematic uncer- 
tainties associated with the prior. We will not consider 
these further. The last parameter, q, arises from the fact 
that the period T is observed but it is the semimajor axis 
that determines the transit probability. In fact, the com- 
panion mass M c has a definite value at the inclination at 
which transits occur, and so this can be determined. We 
explore this latter issue in Section [772] 



7.1. Measurement Uncertainties 

Discussion of transit probabilities — prior or otherwise 
— is often unaccompanied by any discussion of the uncer- 
tainties in the measurements of t he relevant orbital pa- 
rameters. Consider Equation (29), which can be rewrit- 
ten as logP tr = log/ + log(l + r)-ilog(l + (7)-|logp st - 
| logT + logg(e, sin a;) + const. Again, we will ignore the 
uncertainty in / and r here, and we discuss the effect 
of q in the next section. The uncertainty in the tran- 
sit probability ap tr due to uncertainties in the remaining 
quantitities is 



4 (<jt\ 
9 IY/ 



J 9(e 



g(e,smu) 



that is smooth at M c 



M . 



(31) 

where we have assumed that the uncertainties in p*, T, 
and (e, uS) are uncorrelated, which is generally a good as- 
sumption. We note that uncertainty in the transit prob- 
ability depends fairly weakly on the the stellar density 
p*. Furthermore, the period T is usually comparatively 
well-measured, so despite the higher intrinsic sensitivity 
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Period (days) 



Figure 8. Top panels: Posterior probability as determined using the Ida & Lin mass distribution (left) and Mordasini distribution (right) 
for RV-detected planets that have a full list of orbital parameters in the Exoplanets.org Database. The solid grey line shows the simple 
prior transit probability for circular orbits, a star with a solar density, and a planet with a small mass and radius compared to the star. The 
dotted line marks posterior transit probabilities above 10%. Bottom panels: Scale factors for the same planets. The dashed line denotes a 
posterior probability equal to the prior, while the dotted line denotes a posterior that has a 20% boost relative to the prior. In all panels, 
bluer colors mark more massive planets, a plus sign denotes planets whose host stars have R* > 2Rq and X marks planets that have 
Ptr > 0.1 and f a > 1.2; the latter are listed in Table [3] 



Name 



Table 3 

Transiting Planet Candidates 

M c sini(Mj U p) Period (days) M,(Mq) R*{Rq) V (mag) Transit DcptlFl 
1 ~ c 10- 4LJ 



Ptr Ptrfi fa 



61 Vir b 
BD -08 2823 b 
HD 10180 c 
HD 125612 c 
HD 1461 iFJ 
HD 181433~b 
HD 215497 b 
HD 219828 b 
HD 47186 if] 
55 Cnc (FJ 
HD 17156 tF] 
HD 40307 l*f 
HIP 14810 Ff] 
t Boo rj^] 



0.016 
0.046 
0.042 
0.058 
0.024 
0.024 
0.021 
0.062 
0.071 
0.026 
3.30 
0.013 
3.87 
4.17 



4.215 
5.600 
5.760 
4.155 
5.773 
9.374 
3.934 
3.834 
4.085 
0.737 
21.217 
4.312 
6.674 
3.312 



0.942 
0.740 
1.060 
0.902 
1.026 
0.8 
0.872 
1.24 
1.0 
0.905 
1.285 
0.740 
0.990 
1.341 



0.980 
1.279 
1.109 
1.043 
1.130 
1.008 
1.107 
1.468 
1.131 
0.943 
1.507 
0.839 
1.320 
1.418 



1.9 
10 
7.3 
9.0 
6.6 
8.4 
9.1 
8.0 
7.6 
6.0 
8.2 
7.2 
8.5 
1.5 



4.92 x 
8.82 x 
1.05 x 
1.70 x 
5.67 x 
7.03 x 

5.10 x 
9.24 x 
1.79 x 
4.54 x 
5.29 x 
5.33 x 

9.11 x 
8.52 x 



10 
10" 
10 
10 

io- 

10 
10 
10 
10 

io- 4 

10~ 3 
10" 4 
10~ 2 
10" 2 



-3 
-3 



-4 



-3 



0.107 
0.113 
0.103 
0.126 
0.125 
0.119 
0.172 
0.159 
0.130 
0.366 
0.139 
0.109 
0.134 
0.197 



0.081 
0.093 
0.084 
0.105 
0.096 
0.092 
0.131 
0.133 
0.108 
0.284 
0.112 
0.083 
0.099 
0.140 



1.317 
1.216 
1.221 
1.206 
1.302 
1.304 
1.308 
1.203 
1.204 
1.289 
1.240 
1.315 
1.353 
1.405 



^Transit depth (R c / R*) 2 , assuming the minimum mass and a mass-radius relation R c = Rg^Mc/M^) 53 I Weiss et al.|2013 l 
^Rivera et al.|(|2010^ did not definitively rule out transit. 

Candidate using posterior from Mordasini mass distribution. 
d Known transiting planet. 
c Transits excluded. 
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of P tr to it, it generally contributes negligibly to the total 
uncertainty. 

Turning to e and w, we note that the expression for 
c ff ( ej sinu) depends on the covariance of e and uj. For 
simplicity, we consider only the limit in which these pa- 
rameters are uncorrelated: 



a 



where we have dehncd 
4e 2 



E 



Eat + Via 2 , 



■ 2 

sin uj 



(1 - e 2 ) 2 (1 + esinw) 2 



and 



n 



e cos uj 



1 



(32) 



(33) 



(34) 



Figure [9] shows the prefactors E and ft as a function of 
e and uj. 

It is instructive to consider the behavior of E and fl 
is the limit of e <C 1. To first order in e, we have that 
E ~ sin 2 uj — 2e sin 3 uj. For uj — and uj = it, i.e. orbits 
such that the semimajor axis is perpendicular to the line 
of sight, the contribution of the eccentricity to the uncer- 
tainty in P tr vanishes. In addition, for e = 0, there is a 
finite contribution to the uncertainty in P tr for all values 
of uj except and 7r, reaching a maximum of E = 1 for 
uj = 7r/2 or 37r/2; these values correspond to orbits such 
that the semimajor axis is along the line of sight. With 
regards to CI, there is no contribution to the uncertainty 
in P tr to first order in e, and we can therefore expect 
the contribution of the uncertainty in e to dominate over 
that of uj. We conclude that for small e, uncertainties in 
e and uj generally do not lead to large uncertainties in 
Ptr, with the contribution to the fractional uncertainty 
in P tr at most of order a e . 

On the other hand, for orbits with large eccentricity, 
the uncertainties in e and uj can have have a substan- 
tial effect on the uncertainty in the transit probability. 
First consider the full expression for E. It is clear that 
since both terms are positive and the first term in E 
does not depend on uj, E is minimized at uj — for 
all e. Furthermore, E diverges as e —¥ 1 for all uj. 
Next, consider the full expression for ft. The contri- 
bution of uj to the uncertainty vanishes for uj = tt/2 
or 37r/2, regardless of eccentricity. f2 is maximized at 
uj = it + arcsine and uj = 2ir — arcsine, with maximum 
value of fl m ax — e 2 /(l — e 2 ) that diverges as e — » 1. 

7.2. True versus Minimum Semimajor Axis 

In order to estimate the transit probability, one needs 
an estimate of the true semimajor axis. However, one 
measures only the period T and the mass function M.. 
These, together with an estimate of M„, only allow one 
to measure a minimum semimajor axis a m i n . This mini- 
mum semimajor axis can be determined by first using the 
expression for the mass function (Equation 3| to deter- 
mine the minimum companion mass M c ^ m i n for sini = 1, 
and then using this to solve for a using Newton's version 
of Kepler's Third Law (Equation |4| . Since a m i„ is the 
minimum semimajor axis, the transit probability esti- 
mated in this way will be an overestimate of the true 
transit probability. 



Since the expression for a m i„ is complicated, we will 
consider a somewhat simpler estimate for a, 



a 



( GM* 
{ 4tt 2 



\ 1/3 

J T 2 / 3 = a(l + q)- 1 ^. 



(35) 



The relation between the transit probability Pt r ,a esti- 
mated using a,Q and the true transit probability is 



Ptr = Ptr,a„(l + q)^^ ^ P t r,a ( 1 ~ \ 



(36) 



where the rightmost expression holds for q <C 1. Thus, 
for small q, the fractional amount by which one overes- 
timates the transit probability using ao is ~ q/3. For 
a Jupter-mass companion, this is only ~ 0.03% — es- 
sentially negligible. For a companion at the hydrogen- 
burning limit, the error is 2.5%. Of course, one does not 
know the mass of the companion a priori. 

Although the amount by which one overestimates the 
transit probability is generally negligible, the true tran- 
sit probability can nevertheless be determined with no 
approximation. This is accomplished by solving for the 
value of q at the critical inclination for a transit i m in t and 
then using this value of q to determine the true semima- 
jor axis. We proceed by defining 



Q = 



(1 



\2/3 



(1 - cos 2 



) l/2 



(37) 



which we call the "mass ratio function" in analogy to 
the mass function. The true transit probability at the 
critical inclincation is simply the cosine of the critical 
inclination, i.e., P tr = cosi min = Ptr,a {^ + q)~ 1/3 - We 
can then insert this expression for cos i m in into Equation 
(37) and, by defining y = 1 + q, derive a polynomial 



equation in y when i 



(y-i)V /3 



r>2 \ 

tr,a,Q J 



0. 



(38) 



Recall that Q and Ptr,a are observables. This expression 
can be so lved for y with standard techniques. 

shows Ptr/Ptr,a = y^ 1 ^ 3 as functions of 
For typical prior transit probabilities, Pq 



Figure 
Q and P t 



10 



tr, do ■ 



overestimates the actual prior transit probability by less 
than 10% over the entire range of Q. Note that, in the 
limit that q <^ 1, 



q — <3(i - Pt, 



-1/2 



(39) 



In the right panel of Figure 10 for Pt r ,a *C 1, 

Ptr/Ptr,a — 1 _ Q/3. 

8. SUMMARY AND CONCLUSION 

We have derived a general expression for the poste- 
rior probability that a radial velocity companion with 
minimum mass Mo = M c sin i also transits its par- 
ent star. The posterior transit probability depends 
on both the prior distribution of the orbital inclina- 
tion angle i and the prior distribution of the true com- 
panion mass M c . We evaluated this expression for 
power-law distributions of M c and integer power-law 
indices —3 < a < 3, deriving exact analytic expres- 
sions for the posterior transit probility for these cases. 
We found that, for power-law distributions in general, 
the posterior transit probability is well-approximated 



1G 




Figure 9. Contributions to the transit probability error from f2 (top) and E (bottom) as functions of u (left) and e (right) in the 
uncorrelated limit. 



as a scalar multiple f a of the prior probability, with 
specific values /_ 3 = 3/2, /_ 2 = 4/tt, /_i = l,/ = 
2/7r,/i = l/arctanh(cosz min ), / 2 = tanz mi „, and / 3 = 
2/[arctanh(cosz mi „) + cot(i mi „) csc(i m m)], where i min 
is the minimum inclination angle corresponding to the 
maximum companaion mass such that M c>max / 'Mq = 
l/sini mm . For a < 0, f a ~ (-2a/7r) 1 / 2 . 

We then applied our findings to four synthetic but 
physically motivated companion mass distributions. In 
all four cases, we find that a ~ —1 for Jupiters (lOOM^ — 
1O 3 M0), and therefore the prior and posterior probabil- 
ities are very similar. On the other hand, we find a ~ 
-1.5 for Earths and Super-Earths (0.1M© - 10M e ), and 
a ~ -2 to - 2.5 for Super- Jupiters (10 3 M e - 13M Jup ). 
The posterior transit probability for RV planets with 
masses in these regimes is therefore boosted relative to 
the prior, so we may expect more transiting planets in 
these mass regimes. With transit surveys pushing well 
into the Super-Earth regime, this result is encouraging. 

For brown dwarfs (13Mj up — 0.07M Q ), the scale fac- 
tor of the posterior transit probability relative to the 
prior depends in detail on the aridity of the brown dwarf 
desert, as well as on the precise value of the minimum 
mass that is measuremed. In general, the scale factor 
drops as one moves from the Super- Jupiter up through 
the brown dwarf mass regime. In particular, objects 
with minimum mass just above the driest part of the 
brown dwarf desert may have strongly suppressed tran- 
sit probabilities, since such an object is more likely a 



stellar-mass companion being viewed pole-on than a true 
brown dwarf. For stellar companions, (O.O7M — 1M ), 
a ~ 0, and thus the posterior transit probablity is gen- 
erally smaller than the prior probability in this regime. 

Using these mass distributions and the correspond- 
ing scale factors, we estimated posterior transit proba- 
bilities for RV-discovered plane ts in the Exoplanets.org 



Database (Wright et al. 20111. Selecting planets for 
which the posterior transit probability is > 10% and 
> 20% larger than the prior probability, excluding com- 
panions around giants, known transiting planets, and 
planets for which transits have been ruled out, we found 
nine particularly promising transiting planet candidates: 
61 Vir b, BD -08 2823 b, HD 10180c, HD 125612 c, HD 
1461 b, HD 181433 b, HD 215497 b, HD 219828 b, and 
HD 47186 b. 

Finally, we discuss two issues that can subtly affect 
the calculation of either the posterior or the prior tran- 
sit probability. We find that uncertainties in the eccen- 
tricity e and argument of periastron uj generally do not 
yield significant uncertainties in the transit probability, 
provided e is small. However, measurement uncertainties 
in e and uj for companions in highly-eccentric orbits can 
have a significant impact on the transit probability, with 
the uncertainty in e being amplified by a factor that di- 
verges as e — > 1 for all uj, and the uncertainty in u> being 
amplified by a factor that also diverges with increasing e 
for the specific case of uj = it + arcsin(e). Additionally, 
we point out that the semimajor axis typically used in 
determining the transit probability is typically smaller 
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Figure 10. Comparison of the actual prior transit probability, Pt r , and the prior transit probability calculated with ao, Ptr,a , plotted 
against the "mass function ratio" Q (left) and Ptr,a Q (right). For reasonable values of Ptr,a (< 0.5), using ao leads to ovcrcstimation of 
the prior by no more than 10% over the entire range of Q. 



than the true semimajor axis, which leads to an overes- 
timate of the true probability. We demonstrate that this 
overestimate is negligible for systems with small transit 
probabilities, but becomes more significant as the transit 
probability aproaches unity. 

The primary difficulty with estimating the posterior 
transit probablity is that the mass distribution of low- 
mass companions is poorly known. Nevertheless, even 
the theoretical distributions we studied here likely pro- 
vide a more accurate estimate than the naive prior tran- 
sit probability. These distributions suggest that RV com- 
panions with minimum masses in the Earth/Super-Earth 
and Super- Jupiter regimes are promising targets for tran- 
sit follow-up. As measurements of the companion mass 
distribution become increasingly robust, we look forward 
to more frequent use of the posterior transit probability 
— accompanied, of course, by an appropriate discussion 
of assumptions, priors, and sources of uncertainty. 
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